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ABSTRACT 

We address the problem of particle acceleration in solar flares by fast modes which may be excited 
during the reconnection and undergo cascade and are subjected to damping. We extend the calcula- 
tions beyond quasilinear approximation and compare the acceleration and scattering by transit time 
damping and gyroresonance interactions. We find that the acceleration is dominated by the so called 
transit time damping mechanism. We estimate the total energy transferred into particles, and show 
that our approach provides sufficiently accurate results We compare this rate with energy loss rate. 
Scattering by fast modes appears to be sufficient to prevent the protons from escaping the system 
during the acceleration. Confinement of electrons, on the other hand, requires the existence of plasma 
waves. Electrons can be accelerated to GeV energies through the process described here for solar flare 
conditions. 

Subject headings: acceleration of particles — MHD — plasmas — Sumflares — turbulence 



1. INTRODUCTION 

The mechanism of energy release and the process of 
its transfer to heating and acceleration of nonthermal 
particles in many magnetized astrophysical plasmas in 
general, and solar flares in particular, are still matter of 
considerable debate. Recent research shows turbulence 
may play an essential role in these processes. In the case 
of solar flares, it is believed that the energy comes from 
release of stored magnetic energy via reconnection (see 
Priest & Forbes 2000). Turbulence is expected to develop 
since both ordinary and magnetic Reynolds number are 
very large. Both observational evidence and theoreti- 
cal arguments suggest most of the energy is dissipated 
through turbulence (Biskamp 2000; Lazarian & Vishniac 
1999, Lazarian et al. 2004; Shay et al. 2004). Recent 
high resolution observations of solar flares by Yohkoh and 
RHESSI satellites have supplied ample evidence that, 
at least from the point of view of particle acceleration, 
plasma turbulence and plasma waves appear to be the 
most promising agent not only for the acceleration mech- 
anism but also the general energizing of flare plasma (see 
e.g. Petrosian 2007 and Petrosian & Liu 2004 and ref- 
erences cited there). This may also be true in other sit- 
uations (Lazarian et al. 2002, Liu, Petrosian & Melia 
2004). A substantial progress in understanding of in- 
compressible (Shebalin, Matthaeus & Montgomery 1983; 
Higdon 1983; Goldreich & Sridhar 1995, hereafter GS95; 
Matthaeus et al 1996, Cho, Lazarian & Vishniac 2003; 
Biskamp 2003) and compressible MHD turbulence (see 
Lithwick & Goldreich 2001; Cho & Lazarian 2002, 2003), 
as well as MHD-turbulence-particle interactions (Chan- 
dran 2000, Yan k. Lazarian 2002, 2004, 2008, henceforth 
YL02, YL04, YL08, respectively) has clarified many as- 
pects of this problem. 

It was demonstrated that the often adopted Alfven 
modes undergo anisotropy cascade mainly in the direc- 
tion perpendicular to the underlying magnetic field Bq 
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with a Kolmogorov spectrum. As a result they are in- 
efficient in acceleration of particles. Slow modes which 
mix with Alfvenic modes are also negligible for the same 
reason. Fast modes, on the other hand, develop on 
their own, as their phase velocity is only marginally af- 
fected by mixing motions induced by Alfven modes. Ac- 
cording to Cho & Lazarian (2002) fast modes follow an 
isotropic "acoustic" cascade where the velocity scales as 
Vk oc /c" 1 / 4 and in each wave- wave collision a small frac- 
tion of energy equal to v p h/vk is transferred to smaller 
scales. In low [3 medium, the three dimensional energy 
spectrum is (YL02) 

w(k) = nm £ v \-i L -i (WfeM (1) 

with a cascade time scale of 

Teas = (vph/vk^kvk)- 1 = {L/5V)M A \kL)- 1 ' 2 . (2) 

Here k% t j refers to the x,y components of the wave vector 
k, n is the density of the plasma with ions of mass rrii 
(~ proton mass), 5V is the initial perturbation injected 
at the outer scale L = fc m ; n ), v p h = oj/k is the phase 
speed of fast mode with frequency uj and wave vector k, 
and Ma — SV/va is the Alfvenic Mach number for the 
Alfven speed ~ va = B / ^/Airnmi. 

On small scales, the spectrum of turbulence is affected 
by damping. Damping becomes important at a wave 
vector k c when the damping time Tj" 1 becomes compa- 
rable or shorter than the cascading time r cas . Beyond 
this wavevector the turbulence spectrum falls off rapidly. 
In fully ionized plasma, there are basically two kinds 
of damping. Collisional damping is important on scales 
greater than the Coulomb collision mean free path Ac ou i ■ 

4 We discuss in §5 the limitations of the model of turbulence 
in Cho &; Lazarian (2002, 2003). On the basis of weak-turbulence 
theory, Chandran (2005) has argued that high-frequency fast waves 
with fey >> k± interacting with other fast waves generate high- 
frequency Alfven- waves with feii >> kj_. We expect that the scat- 
tering by thus generated Alfven modes to be similar to scattering 
by fast modes that created them. Therefore, within the simplified 
approach adopted in the paper, we do not consider this type of 
interactions. 
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TABLE 1 
Notations in this paper. 



TABLE 2 

The physical parameters of solar flares we adopted. 



a power law index of nonthermal particle distribution 

A(E) acceleration rate 

Eo lower energy limit of nonthermal particles 

k wave vector 

u wave frequency 

k c turbulence cutoff wavenumber due to damping 

L energy injection scale 

8V injection speed 

Ma Alfvenic Mach number SV/va 

n number density of the corona 

N(E) number density of nonthermal particles per energy bin 

No N(E) at the lower energy limit Eq 

T corona temperature 

va Alfven speed 

v particle speed 

fjj cosine of particle pitch angle 

W(k) turbulence spectral energy density 

a ratio of electron mass m e to ion mass 

B magnetic field 

(3 P ratio of gas pressure to magnetic pressure 

P v/c 

Pa v a /c 

r wave damping rate 

Teas cascading time of the turbulence 

T acc acceleration time scale 

T loss energy loss time scale of particles 

T esp escaping time of CRs from the system 

77 cos 9 

Vc va/v 

\coul Coulomb collisional mean free path 

Ay parallel mean free path of CRs 

6 pitch angle of a fast modes 

89 variation of 9 in turbulence 



In solar flares Ac c 



5 x 10 cm I 
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and the relevant scales are shorter so that collisionless 
damping is dominant. Taking into account interactions 
with thermal and nonthermal particles, Petrosian, Yan 
& Lazarian (2006, henceforth PYL06) studied damping 
of fast modes in solar corona condition and showed that 
most the damping is also highly anisotropic. For most 
angles of propagation the inertial range is truncated at 
large scales. Only quasi-parallel and quasi-perpendicular 
waves reach short scales comparable to ion gyroradius 5 . 
For plasma f3 p = P ga s/Pmag ;$ 0.1, the damping due to 
electrons dominates as it's easier for electrons to catch 
up with the waves. The corresponding damping rate is 
(Ginzburg 1961, YL02) 



y/iraPp sin 2 

— — LO- 



COS 6 



exp 



(3 P cos 2 



(3) 



where a — m e /rriH, is the angle between the wave vec- 
tor and the magnetic field. By equating the above equa- 
tion and eq.©, we obtain the cutoff scale of turbulence 
owing to the collisionless damping, 

k - L = - - -, ex P I a —*7i ) ' ( 4 ) 



iraftp sin 9 



(3 P cos 2 



In what follows we shall use this spectrum of fast modes 
to calculate the acceleration and confinement of particles 
by fast modes in solar flare conditions. 

5 The quasi-perpendicular most likely will be damped because 
of magnetic field wanderings (see PYL06). 
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In view of various difficulties with quasilinear theory 
(QLT), a number of nonlinear theories (NLT) have been 
proposed (see Dupree 1966; Volk 1973; Jones, Kaiser 
& Birmingham 1973; Goldstein 1976; Felice & Kulsrud 
2001; Matthaeus et al. 2003; Shalchi 2005). Based on 
particle trapping due to large scale magnetic perturba- 
tions (Volk 1975), we developed a nonlinear formalism 
in YL08 to treat cosmic ray scattering in MHD turbu- 
lence. In view of these progresses, we believe the time 
is ripe to investigate the stochastic acceleration by the 
tested model of MHD turbulence in solar flares. In §2 
wo describe how the nonlinear effects are formulated. In 
§3 and §4 we present results on acceleration and confine- 
ment of the particles, and in §5 and §6 we present a brief 
discussion and summary of the results. Some mathemat- 
ical details are given in the appendix. 

2. NONLINEAR THEORY FOR PARTICLE ACCELERATION 

The interaction of charged particle with MHD turbu- 
lence has been mostly described by the quasilinear theory 
(QLT). While the QLT allows easily to treat the CR dy- 
namics in a local magnetic field system of reference, a 
key assumption in the QLT, that the particle's orbit is 
unperturbed, makes one wonder about the limitations of 
the approximation. While QLT provides simple physical 
insights into particle interaction with turbulence, it is 
known to be an approximation. Recently there has been 
a surge of interest in extending the treatment of this 
problem beyond the QLT. Examples of this are recently 
developed nonlinear guiding center theory by Matthaeus 
et al. (2003), weakly nonlinear theory by Shalchi et al. 
(2004) and second-order quasilinear theory by Shalchi 
(2005a). Most of the analysis so far are confined to slab 
and 2D model of MHD turbulence. In YL08, we extended 
the nonlinear treatment to models of large scale turbu- 
lence and obtained reasonable results for the diffusion 
coefficients. Here, we adopt the same approach for treat- 
ing the particle acceleration for conditions appropriate 
for the solar corona. 

In MHD turbulence, there are basically two types of 
interactions: gyroresonance and transit time acceleration 
(TTD). The resonant condition, for a particle of velocity 
v, pitch angle cosine \x and Lorentz factor 7, is to— k^vji ~ 
nfl/"/ (n — 0, ±1, 2...), where to is the wave frequency, fcy 
is the projection of wave vector along the local magnetic 
field B and f2 = eB/(mc) is the gyro- frequency of the 
particle (charge e and mass m). 

TTD formally corresponds to n = and it requires 
compressible perturbations. When particles are trapped 
by moving in the same speed with waves, an apprecia- 
ble amount of interactions can occur between waves and 
particles. Since head-on collisions are more frequent than 
that trailing collisions, particles gain energies. For small 
amplitude waves, particles should have parallel speed 
comparable to wave phase speed (y p h ~ to/k) to be 
trapped in the moving mirrors. This gives rise to the 
above Cherenkov condition. In general, the momentum 
diffusion coefficient due to interactions (including both 
gyroresonance and TTD) with compressible modes can 



3 



be written as (Yan & Lazarian 2004) 
D pp = nVlm v A (l — [i 2 ) 



d A k 



U B 



Rn(k\\v\\ -uj±nn)J r ? (w),(5) 



where u — k±v±/£l and Ub — B 2 /8tt is the energy den- 
sity of mean magnetic field, Wj_(k) is the kinetic energy 
of the turbulence motions perpendicular to the magnetic 
field. For fast modes in low (3 medium, the velocity per- 
turbations are perpendicular to the magnetic field (Yan, 
Lazarian & Draine 2004), thus Wj_(k) = W(k) (see 
EqQ]). In quasi- linear approximation, particle's orbit is 
unperturbed and therefore the above resonance condition 
should be strictly observed so that Rn(k\\v\\ —u>±ntt) = 
5(h\Vfj, — oj ± nfl). 

However, the assumption of unperturbed orbit is not 
exact, as shown in YL08, the large scale motion (partic- 
ularly compression) changes particles' pitch angle due to 
the conservation of adiabatic invariant. Because of the 
perturbation (especially compression) of magnetic field 
along the particle's trajectory, pitch angle is changing 
according to Av\\/vx ~ (< SB 2 > /B 2 ) 1 / 4 , where v\\ 
and Vj_ are the particle speed along and perpendicular 
to the magnetic field, and SBn is the large scale pertur- 
bation of magnetic strength along the mean field. In this 
case, the resonance function will be broadened (YL08): 



R n (ku vh —Lu±nQ) 



I fcn Attn 



exp 



(kiiVfi — to ± ntt/^y 

kJAV 2 



(6) 

This replaces the 8 function in the QLT formula both for 
TTD (n = 0) and gyroresonance interactions (n = ±1). 
This modification is particularly important for TTD be- 
cause of the existence of a critical pitch angle cosine 
~ va/v given by the Cherenkov resonance condition. 
Physically, TTD take places on all scales while gyroreso- 
nance happens on a local scale. Therefore the large scale 
trapping is more influential to TTD. 

From this we can the obtain the systematic accelera- 
tion rate (see Petrosian & Liu 2004) 

d[vp 2 D(p)} 



A{E) 



Ap 2 dp 



where for isotropic wave distribution, 



1 f 1 

Dip) = ^ ] D p 



(7) 



(8) 



is pitch-angle averaged coefficient. The acceleration time 
is then given by 

r acc = E/A(E) (9) 

3. PARTICLE ACCELERATION BY FAST MODES 

Gyroresonance requires electromagnetic perturbations 
at a particular scale, e.g. k res ~ ~ r ~ . Given 

the parameters we adopt here, k res L > 10 9 for non- 
relativistic electrons. This is certainly beyond MHD 
regime Indeed nonthermal protons can be accelerated 
through gyroresonance with the fast modes. Inserting 
the spectrum of fast modes from equations (fTJ) into equa- 
tion©, we get 



D, 



G _ vV^M 2 A m 2 v 2 A (l - M 2 ) 
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exp 



( M _24±l ) 
' v xr/H ' 



(10) 



Nevertheless, our calculations shows the overall contribu- 
tion from gyroresonance is much smaller than that from 
TTD. In fact for the adopted parameters, there is no 
acceleration from gyroresonance for protons with ener- 
gies less then <4MeV because of the damping of fast 
modes. For higher energy protons, acceleration rate from 
gyroresonance is 2-3 orders of magnitude lower than that 
from TTD. 

On the other hand, TTD interactions can happen on 
all scales or wavevectors k and appears to be the dom- 
inant acceleration mechanism for protons and the only 
channel to direct the energy from MHD turbulence to 
electrons. The corresponding momentum diffusion, in 
the NLT limit, can be obtained in a similar way, 
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where Afj, = Avu/v, R = vkj/Q is the dimensionless 
rigidity, M A is the Alfvenic Mach number, and 77 = cos 9. 
Combining equation (TTTj) with equation (|7I8|) we get the 
acceleration rate 
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Comparison of NLT with QLT: The QLT momentum 
diffusion and direct acceleration rates due to TTD are 
obtained by adopting the delta function as the resonance 
function R n (kuvu — lo) — b~{k\\y\\ ~ UJ ) : 



nM\rn 2 v 2 A v 2 . . . , 
D m = A (1 - M ) / V d V 
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For electrons and low energy ions (< lOOKeV), in mag- 
netic field of B — 10 2 G, vj_/fl e < c/Q e ~ 14cm, so 
for most angles k±v±/Q << 1 and the Bessel func- 
tion can be replaced by the first order approximation 
J n {x) ~ (x/2) n /n\. Thus inserting k c given by eq.((4]) 
into equation (fTI)) . we obtain for the QLT limit 
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Fig. 1. — The acceleration tim es based on NLT (dashdot line, 
Eg 112(1 and QLT (solid line, Eg |14j l, The lower thick lines represent 
the result for electrons and the upper thin lines are for protons. 
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where $/ is the imaginary error function. The NLT and 
QLT result are compared in FigHJ We see that the two 
results qualitatively agree with each other so that for 
most cases one can use the acceleration rate using QLT. 

Effects of field line wandering: Unless the medium is 
strongly magnetized {(3 P < 0.1), field line wandering can 
not be neglected (YL04). As shown in PYL06, field line 
wander about < 15° around 90°. Since for high velocity 
particles, rj w va/(v(j,) -C 1 and the quasi-perpendicular 
modes dominate. Thus this angle variation around 90° 
must be taken into account. The damping rate should 
be averaged over the range 90° — 59 ~ 90°. Equating 
cascading rate (see Eqj2]) and the damping rate (Eql3]) 
averaged over 59 (PYL06) 
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(16) 



and combining it with the variation of angle 9 due to field 
line wandering 59 ~ (M^/fcL) 1 / 6 , one gets the averaged 
damping wavenumber k c . For Ma = 0.3, it has only so- 
lution for f3 p > 0.04. For (3 p = 0.1, k c L ~ 16. In FigEJ 
we compare the acceleration rate with and without field 
line wandering. As can be seen, the acceleration is de- 
creased in the case of (3 P > 0.04. This can be explained 
as follows. According to the Landau resonant condition 
fciiun ~ u>, the higher the energy, the closer to 90° the 
resonant wave vector. As shown in Eq. (|16p . the rapid 
increase of the cutoff wavenumber near 90° is smeared 
out due to the field line wandering. Acceleration rate is 
accordingly modulated for particles with different ener- 
gies. This result shows that field line wandering has a 
noticeable impact on the acceleration and must be ac- 
counted for. Because of the field line wandering, accel- 
eration is substantially reduced when (3 P > 0.04. For 
j3 p < 0.01, effects of field line wandering is negligible due 
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Fig. 2. — The influence of field line wandering to TTD acceler- 
ation by fast modes. The solid and dotted lines are the averaged 
accelerate times taking into account field line wandering; the dot- 
ted lines are the results without accounting for field line wandering. 
The acceleration rate is decreased due to averaging because field 
line, wandering increases damping, at large angles. „ . „ , , . 

to the dominance of magnetic energy. Tn figures 13141 we 
compare the acceleration time due to TTD interaction 
with other times for electrons and protons. We see an 
increase of acceleration rate at ~ lOKeV for electrons 
and ~ lOMeV for protons because of field line wander- 
ing. The acceleration times drop suddenly because the 
wave damping drops dramatically when the correspond- 
ing 9 approaches 90° . The fact that the acceleration rate 
increases substantially at (3 P < 0.04 sets a critical value 
for the plasma (3 p < (3 cr provided that the stochastic 
acceleration by MHD turbulence is the dominant mech- 
anism for generating high energy particles. Here [3 cr is 
approximately determined by Ma through the Eq. (|16[) . 

4. CONFINEMENT OF PARTICLES 

For efficient acceleration particles should be confined 
within the system for a sufficient duration. This is par- 
ticulary nevessary for TTD where particles gain energy 
only in the parallel direction. They need to be isotropized 
in order to reach high energies. The isotropization is 
dominated by Coulomb collisions for low energy protons 
(< 100 kcV). The Coulomb scattering rate drops quickly 
with the increase of energy. For higher energy particles, 
the scattering is owing to the interaction with turbu- 
lence. The same turbulence that accelerate the particles 
statistically also cause pitch angle scattering and thereby 
provide confinement for these particles. 

We showed in YL08 that quasilinear approximation is 
not valid for scattering from TTD and resonance broad- 
ening due to large scale trapping must be taken into ac- 
count. This removes the divergence at the 90° and leads 
to a finite mean free path, 
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where the pitch angle coefficient due to TTD is 
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electron acceleration by TTD ((3=0.01) 
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Fig. 3. — times involved in the acceleration of electrons. The solid 
line refers to the acceleration time. The dashed line represents en- 
ergy loss time due to Coulomb collisions and synchrotron radiation. 
The dashdot lines give the escaping times. Scattering with MHD 
modes (lower curve) is ineffective; the actual confinement is dom- 
inated by the interaction with the quasi-parallel whistler modes 
(upper curve). The cross of accelerate time and energy loss time 
determines the threshold at the lower energy end. The acceleration 
continues till the cross at the higher energy end. 
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Fig. 4. — Same as Fig(3] but for protons. The acceleration 
starts with super-Alfvenic particles (~ 2(T/1Q 7 K )//3 p KeV). The 
confinement is owing to the scattering by the MHD modes (dashdot 
line). The acceleration continues to O.lGcV, beyond which the 
escaping becomes much faster than the acceleration. 
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TTD dominates the scattering for most pitch angles. 
For high speed particles with small pitch angles, how- 
ever, TTD interaction is negligible as indicated from the 
Chenrenkov condition and gyroresonance takes over. We 
showed in YL08 that QLT provides a good approxima- 
tion for gyroresonance since it is a local interaction on 
small scales. We adopt here the QLT result for the pitch 
angle scattering due to gyroresonance, 

D«= 27 r0 2 (l-/, 2 ) £ d?k^> 

a(jfei|«n -u±si) (v + ^) 2 A {uf ■ (19) 

The confinement time of the particles can be approxi- 
mated by (see Petrosian & liu 2004) 
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This time scale is also shown in Figs l3l4l The thresh- 
old for TTD interaction is vu > v A as indicated by 
the Cherenkov resonance condition. That's why the 
confinement time for the protons L 2 /(Xv) rises around 
Ek — lOOKeV and drops on high energy end (see FigH]). 
Acceleration continues up to ~ lOMeV, beyond which 
the confinement time becomes much smaller than the 
acceleration time. 

Indeed pitch angle varies with the variation of magnetic 
field due to conservation of adiabatic invariant. The large 
scale slow variation of magnetic field (compared to the 
gyration of particles), however, is inefficient in accelerat- 
ing the particles (see Cho & Lazarian 2006). Therefore, 
while for acceleration QLT and NLT estimates are com- 
parable, scattering and spatial diffusion has to be studied 
with nonlinear approach to remove the 90° degeneracy 
and to obtain finite mean free paths (see YL08). 

For electrons, lack of gyroresonance (except for high 
energy electron with r res > A:" 1 ), the interaction with 
MHD modes is not enough to confine the particles. A 
cutoff of the spectrum > lOOkeV would occur if there are 
no other interactions. Higher frequency whistler modes 
have been considered as one candidate to interact with 
electrons (Wentzel 1976; Melrose & Brown 1976; Be- 
spalov et al. 1991; Stepanov et al. 2007). However, 
due to their anisotropy (Cho & Lazarian 2004), the in- 
teraction with whistler turbulence is very inefficient sim- 
ilar to the case of CR scattering by Alfvenic turbulence 
(Chandran 2000; YL02). On the other hand, parallel 
propagating whistler waves, can be generated by kinetic 
instabilities, e.g., beaming instability, anisotropy insta- 
bility, cyclotron instability (Akhiezer et al. 1975; Tsy- 
tovich 1977; Neubert & Banks 1992) and be an efficient 
agent to scatter electrons. Through the TTD accelera- 
tion, an anisotropic distribution of particle with respect 
to the magnetic field can be generated, which can induce 
cyclotron instability (Lazarian & Beresnyak 2006). Leak- 
age during solar flares is also possible to create beaming 
instability. Scattering by the whistler modes limits the 
bulk speed of the particles close to the phase speed of 
the whistler waves. The waves are generated under the 
cyclotron resonance u = kv\\ — il e /~/ e . Combined with 
the dispersion relation of whistler modes 

" = F + Ke/C) 2 ' ^ 

from which we can get the resonant wave number 
kresile)- The group velocity of the resonant modes 
is then obtained by inserting it into v g (^ e ) = dw/dk, 
which is a function of the electron energies. The confine- 
ment time will be approximately L/v g . For high energy 
electrons (Ek > lOMeV), the resonant wave modes get 
smaller than the gyrofrequency of thermal ions and the 
group velocity v g ~ va for the wave modes moving par- 
allel to the magnetic field. Fig[3] shows that the parallel 
propagating whistler modes are able to confine the elec- 
trons during the acceleration. 

The whistler waves are subjected to both Coulomb 
damping and collisionless damping. For electrons with 
energies < lGeV, however, whistler wave can be excited 
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and provide effective scattering for the electrons (Dor- 
man 1996). The detailed study is beyond the scope of 
this paper. 

5. DISCUSSION 

In this paper we have discussed the effects that the 
cascade of turbulence has on acceleration and heating of 
Solar corona. There the energy is injected at large scales 
much larger than any plasma scale concerned and this 
justifies a magnetohydrodynamic treatment of the large 
scale motions. Due to recent insights into the physics 
of MHD cascade and its interaction with charged parti- 
cles we reduced a complex problem of acceleration and 
heating to a more manageable problems of interactions of 
Alfven, slow and fast modes with plasma and energetic 
particles. 

Alfvenic turbulence is inefficient in scattering and ac- 
celerating particles because of its anisotropy (Chandran 
2000, YL02). Fast modes, instead, have been identified 
as the MHD turbulence modes dominating the interac- 
tion with cosmic rays (YL02,YL04). In this paper, we 
apply the result to the stochastic acceleration in solar 
flares. 

We assume that the MHD turbulence is strong, i.e. 
that the critical balanced condition is satisfied for 
Alfvenic modes. Therefore, to describe fast mode of 
MHD turbulence we appeal to the results by Cho & 
Lazarian (2002, 2003) on the scaling and coupling of 
Alfvenic and fast modes. The corresponding papers 
claim the isotropy of the fast modes. One may expect to 
see deviations from isotropy, however. For instance, slab 
Alfven modes created by streaming instabilities are sub- 
ject to non-linear damping by the ambient Alfvenic tur- 
bulence (YL02, YL04, Farmer & Goldreich 2004, Lazar- 
ian & Beresnyak 2006). Beresnyak & Lazarian (2008) 
showed that the corresponding damping of the quasi- 
parallel modes depends on the angle between their k 
vector and the direction of magnetic field. One might 
expect quasi-slab fast mode to be subject to a similar 
damping by strong Alfvenic turbulence. This, however, 
has not been demonstrated so far. 

If at the injection scale SB <C B, the Alfvenic tur- 
bulence is weak (see Galtier et al 2000) and develop a 
cascade with fcj| = const. The interaction of the weak 
Alfvenic turbulence with other modes can be very dif- 
ferent from that of the strong Alfvenic turbulence. For 
instance, it was shown by Chandran (2005) that fast 
modes develop anisotropy (i.e. the energy in the quasi- 
slab modes is reduced) owing to their interaction with 
Alfven modes in the weak regime. However, such the 
Alfvenic weak turbulence has a limited inertial range (see 
discussion in Cho, Lazarian & Vishniac 2003) and at suf- 
ficiently large k transfers into a strong turbulence. More- 
over, magnetic reconnection in Solar Flares should pro- 
duce perturbations SB ~ B, which should induce strong 
turbulence from the very beginning. 

Strong Alfvenic turbulence is characterized by k± ^> 
fey with the GS95 relation fc|| ~ fcj_ defining a cone in 
the Fourier space where most of the turbulent energy re- 
sides. However, the aforementioned relation should not 
be understood too literally. The energy outside the cone 
is not zero and the modes with k\\ > k± are present. 
Such Alfvenic modes are weakly interacting even being 
a part of the strong Alfvenic turbulence. Therefore, the 



Chandran (2005) model is applicable to them. Neverthe- 
less, according to Cho, Lazarian & Vishniac (2002) the 
energy in these modes is exponentially reduced 6 . There- 
fore, assuming that the Alfvenic turbulence is injected 
at a scale much larger than the typical gyroradius of the 
energetic particles, we did not consider the anisotropies, 
introduced by the process of the fast wave cascading in- 
duced by Alfvenic modes with large fcy and small k± (cf. 
Chandran 2005)." 

The model for MHD turbulence that we adopted in the 
paper is the turbulence where the back-reaction of ener- 
getic particles is limited to changing the cut-off of the 
turbulence. A more fundamental modifications of turbu- 
lence are conceivable, however. For instance, Lazarian & 
Beresnyak (2006) argued that compressions of the fluid 
of energetic particles may result in the gyroresonance in- 
stabilities that can induce an additional quasi-parallel 
component of the Alfvenic waves. If true, these waves 
would interact with fast modes in the manner described 
by Chandran (2005), which would affect the fast mode 
isotropy. However, if substantial portion of energy re- 
sides within these quasi-slab modes, their major effect 
will be the direct gyroresonance acceleration of energetic 
particles. We felt that the modification of MHD turbu- 
lence arising from the instabilities within energetic parti- 
cles, e.g. by the process in Lazarian & Beresnyak (2006), 
is beyond the scope of the present paper. 

Apart from the issue of isotropy, there are potential 
issues related to the exact scaling of fast modes. One 
dimensional numerical simulations in Suzuki, Lazarian & 
Beresnyak (2007) indicate that fast modes may develop 
a shock-like cascade, which differs from the finding in 
3D MHD calculations in Cho & Lazarian (2002, 2003). 
We adopted the model from the latter works, but wait 
for higher resolution 3D numerical runs to rectify the 
scaling. 

In addition to being strong, MHD turbulence that we 
considered was balanced in the sense that the equal flux 
of energy was assumed in every possible direction. The 
properties of imbalanced turbulence (see Cho, Lazarian 
& Vishniac 2002; Lithwick, Goldreich & Sridhar 2007; 
Beresnyak & Lazarian 2007, Chandran 2008) can be very 
different from the balanced one. However, we expect the 
flow of Alfven waves, which constitute the weak turbu- 
lence, to be subjected to reflection within a Solar corona 
environment that we deal with. As the result we expect 
only marginal imbalance for the problem that we deal 
with. 

A threshold for the TTD acceleration is v > va set 
by the Cherenkov resonance condition. In the low (3 P 
environment, the thermal protons can not be accelerated. 
The low energy threshold would be ~ 2(T/10 7 )//3 p KeV. 
Thermal electrons, instead, can have TTD interactions 
unless the plasma beta is too low (3 p < m e /m p . Our 
calculations show that TTD acceleration dominates over 
gyroresonance for the energy range we consider. The 
acceleration efficiency decreases with the plasma (3 P as 
damping of the fast modes increase with f3 p . 

Particle acceleration rate depends on the wave spec- 

6 Because of this exponential reduction the rates of Alfvenic 
scattering in YL02, which was the first study to take into account 
the modes outside the GS95 cone, were still grossly sub-dominant 
to the fast modes. 
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trum and the wave damping rate are partially determined 
by the particle spectrum. In general, it is required that 
a self-consistent treatment of the evolution of turbulence 
and particle acceleration. The calculations above as- 
sumes that the spectrum of the turbulence is determined 
by the cascade and damping by thermal particles. Our 
numerical calculation of these integrals shows that up to 
~ 10% of total energy of turbulence is being transferred 
to nonthermal particles for the given set of parameters. 
In most cases, the back reaction of nonthermal parti- 
cles is negligible as the damping due to the interaction 
with these nonthermal particles is smaller than thermal 
damping taking into account field line wandering. In 
some large flares, however, a large amount of particles 
need to be accelerated from the thermal reservoir to en- 
ergies 3> k B T (PYL06). In this case, the damping by 
nonthermal particles can be significant and one needs to 
solve the coupled equation of the evolution of turbulence 
and particles simultaneously. 

Isotropization is important for the process of accelera- 
tion. If there is not enough isotropization, the accelera- 
tion by TTD will stop quickly as only parallel velocity is 
increased during the process. Moreover, without enough 
scattering, particles will leave the system before they get 
accelerated. Scattering by fast modes is shown to be 
adequate for isotropization of protons. Different from 
the acceleration, the scattering by gyroresonance, even 
smaller than the TTD scattering rate, plays an essential 
role in determining the scattering of parallel moving par- 
ticles, which can be a substantial portion of the particles 
because of the TTD acceleration. 

For electrons, due to their small gyro-radii, gyroreso- 
nance does not occur with MHD turbulence except for 
those high energy electrons. The actual scattering can 
happen with the plasma turbulence. The work by Cho & 
Lazarian (2004) shows that whistler modes are even more 
anisotropic than Alfven modes. We know gyroresonance 
is very inefficient with anisotropic turbulence. Analo- 
gous to the MHD regime, parallel propagating modes 
may be generated by kinetic instabilities (see Tsytovich 
1977) and be a candidate for interaction with electrons. 
The TTD interaction itself may induce gyroresonance in- 
stability through creating an anisotropic distribution of 
particles with respect to the magnetic field (Lazarian & 
Beresnyak 2006). Our estimate shows that the whistler 
waves (for moderate energy) and Alfven wave (for high 
energy) generated by the anisotropy instability can pro- 
vide effective scattering and isotropization for the elec- 
tron acceleration. 

Acceleration of particles by fast modes were previ- 
ously studied by a number of authors, including Miller, 
Larosa & Moore (1996), Schlickeiser & Miller (1998). In 
their studies, turbulence is assumed isotropic with cither 
Kolmogorov or Kraichnan spectrum. Although coupled 
equations of wave and particle evolutions were solved in 



Miller, Larosa & Moore (1996), we feel that the one di- 
mensional treatment they adopt is problematic, as the 
damping caused by the TTD interaction with particles 
is anisotropic. In this paper, we start with a more phys- 
ically motivated and numerically tested of turbulence. 
We deal with fast modes, which are subject to much 
more efficient linear dissipation. On the small scales, 
however, because of the dissipation above, fast modes 
develop anisotropy. In our treatment, this anisotropy 
is strongly affected by field line wandering which is de- 
termined by the Alfvenic Mach number and plasma (3 
and the efficiency of the acceleration is substantially in- 
fluenced accordingly. In addition, scattering and con- 
finement are treated using the nonlinear theory we have 
recently developed (Yan & Lazarian 2008). We believe 
that the approach that we developed is applicable be- 
yond pure Solar Flare problems, e.g. it can be modified 
to study turbulent acceleration in the medium within 
clusters of galaxies (see Brunetti & Lazarian 2007). 

6. SUMMARY 

We discuss stochastic acceleration in solar flares. The 
dominance of thermal damping makes it possible to de- 
couple acceleration and turbulence evolution. This sim- 
plifies the problem substantially so that we don't have 
to rely on simulation to resolve the problem. Our study 
can be summarized as below: 

• Fast modes are sufficient to accelerate super- 
Alfvenic particles to high energies. Electrons can 
reach GeV energies through the process. 

• Both nonlinear and quasilinear approaches show 
TTD dominates the acceleration by MHD modes. 
For the acceleration, QLT result qualitatively 
agrees with the nonlinear result and appears to be 
a good approximation. 

• Confinement of protons can be realized through 
both gyroresonance and TTD with fast modes. 
Electrons, however, needs plasma perturbations to 
be confined. 

• Magnetic field line wandering reduces the acceler- 
ation efficiency in weakly ionized plasma and this 
sets an upper limit of plasma [3 P for solar flares. 
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APPENDIX 

ENERGY LOSS RATE 
In the paper, the acceleration time r acc = E/A(E) is compared with the energy loss time: 



T loss - E/E loss - 4r 2 (7meCln A/f3 + i?o 2 /3 2 7 2 /9/meC) ( A1 ) 



where tq = 2.8 x 10 13 cm is the classical electron radius and In A = 20 in our regime (see Sturrock 1994). The ion loss 
is mainly due to Coulomb collisions with electrons and protons (Post 1956; Ginzburg & Syrovatskii 1964). The loss 
due to electron-ion collisions is (Petrosian & Liu 2004) 
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At low energies, the loss is dominated by ion-ion collisions, 
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